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OPTICAL  CONSTANTS  OF  GERMANIUM 


ABSTRACT 

An  analysis  has  been  made  of  the  effects  of  change  of  refractive  index  of  ger¬ 
manium  on  radiometric  calibration;  this  analysis  shows  that  the  index  cannot  change 
by  more  than  0.61  if  2%  accuracy  is  to  be  maintained — unless  calibration  proce¬ 
dures  are  used  to  account  for  this  change.  Several  methods  of  measuring  index 
have  been  investigated;  as  a  necessary  preliminary,  the  effects  of  changes  in  tem¬ 
perature  and  pressure  on  index  were  evaluated.  For  the  usual  atmospheric  condi¬ 
tions,  variations  of  50°c  or  so  tore  will  cause  changes  of  no  more  than  ±0.00001 
of  the  refractive  index  of  air.  If  temperature  or  pressure  measurements  can  have 
errors  as  large, as  so°c  or  SO  torr,  then  the  error  or  uncertainty  in  index  can  be 
as  large  as  10*5.  similarly,  the  sample  must  be  measured  to  ±0.04°C. 

Each  of  nine  possible  techniques  has  been  analyzed  in  terms  of  the  precision 
necessary  in  the  measurement  of  each  of  the  independent  variables.  Deviation 
techniques  are  discussed  in  Section  3.1,  reflection  techniques  in  Section  3.2,  inter¬ 
ferometric  methods  in  Section  3.3,  and  a  sphere  method  in  Section  3.4.  The  re¬ 
quired  measurement  accuracies  are  summarized  in  Table  I.  Transmission  meas¬ 
urements  made  on  a  Sample  Submitted  by  the  Sponsor  are  also  described  and 
analyzed. 


1 

INTRODUCTION 

The  nimbus  satellite,  being  developed  by  NASA  and  the  Weather  Bureau  to  determine  basic 
and  dynamic  properties  of  our  atmosphere  and  the  development  of  weather  conditions,  includes 
in  its  instrumentation  a  grating  spectrometer  and  infrared  radiometers  which  use  germanium- 
immersed  thermistor  bolometers.  Both  types  of  instrument  are  meant  to  measure  with  a 
radiometric  accuracy  of  2%  or  better.  The  NIMBUS  satellite  will  undergo  the  usual  temperature 
variations  of  orbiting  vehicles:  a  variation  of  about  -80°C  to  +60°C.  Accordingly,  it  is  of  interest 
to  know  how  the  refractive  index  of  germanium  and  the  absorption  vary  with  temperature. 

.  The  approach  to  solving  this  problem  has  been  to  investigate  the  many  different  methods  of 
measuring  refractive  index  to  determine  which  will  provide  the  most  reliable,  accurate,  and 
economical  answers.  A  comprehensive  literature  search  on  the  properties  of  germanium  was 
also  undertaken,  and  some  measurements  were  made  of  the  transmission  of  poor  germanium 
samples.  This  report  is  divided  into  three  main  categories:  the  evaluation  of  different  methods 
of  measuring  the  refractive  index  of  germanium,  transmission  measurements,  and  the  bibllog- 
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REFRACTIVE  INDEX  MEASUREMENTS 

This  section  omits  a  discussion  of  the  more  usual  Simple  techniques:  Pulfrich,  Abbd, 
Optical  path  measurement,  and  critical  angle  measurement.  The  more  precise  techniques 
have  been  considered  in  greater  detail;  these  include  the  following:  variation  deviation  tech* 
niques,  interferometric  means,  and  reflective  measurements  like  those  of  Robinson  and  others, 
the  analysis  is  based  on  a  desire  to  measure  index  to  o.ooooi;  the  results,  however,  are  in  a 
more  general  form;  one  that  gives  the  required  precision  in  a  measured  value  for  any  accuracy. 

Although  the  extremely  accurate  measurement  of  refractive  index  is  an  important  academic 
problem,  it  is  also  useful  to  investigate  just  how  precise  the  index  measurement  must  be  to 
obtain  a  2%  radiometric  accuracy.  This  analysts  is  also  included  below. 

2.1.  METHOD  OF  ERROR  ANALYSIS 

The  refractive  index  is  a  function  of  several  variables: 

n  »  n(x,  y,  .  . .  ). 

H  one  of  these  variables  changes  by  an  unknown  amount,  it  can  be  written  as  its  original  value 
..  plus  the  change.  Thus  the  index  is  Written  as 

n  +  An  ■  n(x  +  Ax,  y  +  Ay,  .  .  .  ), 

The  change  in  index,  An,  is  written  as 

_  An  «  n(x  +  Ax,  y  +  Ay,  .  .  .  )  -  n(x,  y,  .  .  .). 

Clearly,  this  is  exactly  the  expression  for  the  differential  of  a  function  of  several  variables 
before  the  limits  are  taken.  This  is  the  basis  for  using  differential  calculus  in  error  analysis. 
However,  the  negligibility  of  second-order  difference  terms  must  be  verified  by  one  means  or 
another.  For  example,  the  change  in  total  radiant  omittance  of  a  body  should  be  written 

AW  wq{(*  At)  (T  +  AT)4  -  oeT4 

'  When  expanded,  this  is  written 

AW  ■  o(e  +  At)[4T3  AT  +  OT^AT)2  +  4T(AT)3  ♦  (AT)4]  +  oT4A« 

1  <  '• 

The  terms  which  include' AT  to  a  power  greater  than  one  and  terms  like  acat  must  be  small 
compared  to  AT  and  At  terms.  Then  the  expression  obtained  by  differentiation  is  valid: 

dW  -  «T3AT  ♦  0T4At 
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if  the  second-order  differences  are  small,  calculus  is  valid  and  the  chain  rule  applies: 


.  3n  .  dn  . 
dn'3^dx  +  3?dy  + 


The  relative  error  can  be  written 


dn 

n 


2.2.  EFFECT  OF  INDEX  CHANGE  ON  RADIOMETRIC  ACCURACY 

In  order  to  analyze  the  effect  on  radiometric  accuracy  of  a  change  in  refractive  index,  one 
can  start  with  the  basic  formula  for  ir radiance  (assuming  no  absorption).  Radiance  in  an  optical 
system  is  constant  if  there  are  no  losses.  Assuming  that  the  only  loss  is  at  the  front  surface 
of  the  immersion  lens,  the  irradiance  of  the  detector  is  given  by  H  *  N(1  -  R)u>.  The  change  is 
obtained  by  differentiation: 


t 


dH  =  N[(l  -  R)db>  -  todR] 

d«  dw  dR 
H  o>  *  1  -  R 


Thus  there  are  two  influences:  change  of  reflection  loss,  and  change  of  solid  angle.  The  re¬ 
flectivity  error  of  an  uncoated  sample  is  found  by  assuming  that  the  radiation  is  incident  nor¬ 
mally.  Thus  the  equations  for  reflectivity  and  its  errors  are 


dR  _  4dn(n  -  l)(n  +  1)^ 
’  R  (n  +  l)®2n 


*  2dn 


(n-1) 
n  +  1 


The  error  in  solid  angle  due  to  a  change  in  critical  angle  is  obtained  by  considering  the  change 
in  critical  angle  as  a  function  of  temperature.  The  critical  angle  a  is  obtained  from  Snell's 
.  Law  by  assuming  that  the  refracted  angle  is  n/2: 

sin  a  s  n'/n 

Here  n'  is  the  index  of  selenium  and  n  is  that  of  germanium.  The  error  obtained  by  differentia¬ 
tion  is  found  as  follows: 


cos  a  da 


ndn'  -  n'dn 
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im'  ti\  \  J 
^n  d  J  cos  a 


The  solid  angle  Is  proportional  to  the  square  of  the  critical  angle,  so  that  one  has  the  following 


relation: 


The  relative  error  is  found  by  division: 


o>  =  kor 


dui  =  Skadar 


da)  2  da 


By  Combining  the  above  equations,  the  relative  error  Can  then  be  computed  to  be 


(dn!  n'dn\ 


,dn’  dn 
■  n’  ’  ii 


1/2  arc  Sin  ^  K(n/n') 

Numerical  evaluation  is  now  possible.  If  n  is  4  and  n*  is  2.4  (for  germanium  and  Selenium 
respectively),  then 

m  =  0.4  dn'  -  0.25 dn  +  2dn(n  -  1)  „  „„  dn,  +  „  63  d„ 

H  1/2  arc  sin  0.6  V25/9  *  1  it  +  1 

Therefore,  a  total  change  of  0.02  requires  knowledge  to  about  0,01  in  the  refractive  index  of 
both  the  germanium  and  the  substrate. 


2.3.  EFFECTS  OF  TEMPERATURE  AND  PRESSURE  CHANGES  ON  INDEX  CHANGES 

One  of  the  basic  problems  in  the  analysis  of  methods  for  index  measurement  is  an  evaluation 
of  the  effects  of  the  uncertainty  in  knowledge  of  the  temperature  and  pressure  of  air.  The  equa¬ 
tion  relating  the  index  of  air  to  temperature  and  pressure  [1,  2]  is  given  as 

P(l  +0.P)(1  +  I5a) 

/_  «  \  _  /_  .  I  _ 


("l  *  l)  *  (n15,780  ‘  760(1  +  760d15)(l  +  aT) 
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0t  =  (1.05  -  0.015  T)  10*° 

0lg  =  0.813  X  10*6 
a  =  0.0037 

n, .  ...  =  n'  Is  the  index  of  air  at  15°C  and  760  mm  Hg  . 

lOjioO 

Substitution  yields 

nj  *  1  »  0.001367  (»•  -  1)[P  +  (1.05  -  0.015T)10*6  P2]  (1  +  0.0037T)'1 

The  uncertainty  in  the  index  of  air  with  respect  to  temperature  is  found  by  the  usual  technique 
of  differentiation.  The  partial  derivative  with  respect  to  temperature  IS 

8nl  ,  (n'  *  1)(0.001387){(1  +  0.0037T)(-0.015  x  10'6P2)  -  (0.0037)  [(P  + 1.05  -  0.015T)10'6P2]} 
3T  (1  +  0.0037T)2 

Assuming  as  before  that  T  IS  373°K,  P  is  760  mm  Hg,  and  n'  -  1  is  380  x  10  °,  the  resulting 
-7  - 

value  is  1.936  x  10  .  The  variation  of  index  with  respect  to  pressure  is  found  as  follows: 

W  *  (O'i\10°6S37T~)  t1  +  2(105  -  °-015T)  10'6pJ 

If  it  is  assumed  that  T  *  373°K  and  P  »  760  mm  Hg,  as  before,  the  numerical  result  is  1.62  x 
-7 

10  .  The  total  change  is  given  by  the  total  derivative: 

dnj  •  1.936  x  10_7dT  +  1.62  X  10’7dP 

Pressure  and  temperature  variations  of  less  than  about  50°C  or  50  mm  Hg  will  cause  an  error 
of  no  more  than  ±0.00001. 

A  similar  analysis  can  be  applied  to  the  sample  temperature.  Since  dn-  »  ndn,  +  n.dn, 
n.  dn,  .4  4  4 

where  n  *  — ,  n,  must  be  known  to  2.5  x  10  °.  Since  -rsr  is  about  2.88  x  10  from  earlier 
nl  1  dn.  dT 

measurements  [3]  and  since  -g^r  is  around  three  orders  of  magnitude  smaller  than  this,  then  the 
temperature  accuracy  needed  is  determined  by  the  change  of  index  of  germanium  with  tempera* 
i  ture: 


dn2«  2.68X10  dT 
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So  fig  will  be  accurate  to  ±0.00001  only  within  a  temperature  range  of  about  0.037°C.  Alterna¬ 
tively, 


1  3n2  n2  8ni 

l  nj 


-  2.68  x  10'4dT  «  4  x  1.936  x  l(f7  dT 
*  2.68  x  10"4df 


The  uncertainty  in  temperature  can  be  no  greater  than  1  part  In  27,  or  0.037. 
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DISCUSSION  AND  EVALUATION  OF  VARIOUS  REFRACTIVE  INDEX 
MEASUREMENT  METHODS 


3.1.  DEVIATION  METHODS 

Methods  which  measure  the  deviation  of  a  beam  of  light  from  some  known  direction  are 
among  the  oldest  and  most  accurate  known.  The  classical  method  of  minimum  deviation, 
discussed  in  every  good  text  on  optics,  Is  notable  for  its  simplicity  and  elegance.  It  is  not 
necessary  to  obtain  minimum  deviation  for  the  determination  of  index;  in  some  applications  it 
is  useful  not  to  do  so  (the  technique  then  may  simply  be  called  "prism  deviation").  Also,  some 
simplification  is  achieved  by  measuring  deviations  for  two  prism  positions.  There  is  then  an 
almost  endless  number  of  possible  variations  of  the  techniques;  all  sorts  of  geometric  figures 
may  be  used.  Minimum  deviation,  prism  deviation,  and  two-position  deviation  will  be  discussed 
in  this  section. 

3.1.1.  MINIMUM  DEVIATION  METHOD.  By  the  usual  manipulations  [4]  one  obtains  the 
following  expression  for  minimum  deviation  (Figure  1  defines  the  symbols): 

nj  sin  !/?(®m  *  ®) 

n2  “  sln"<i72 

The  total  derivative  is  given  by 


dn2 


3n2 

8nj  **1  +  86  dflm 


dn2 

*lada 


The  partial  derivatives  are  as  follows: 
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♦ 


3n2  sin  i/2(8m  +  a) 

3n^  3  2  Sin  a/2 

9n2  nj  cos  1/2(9 m  +  a) 

36_  ”  2  sin  cr/2 

m 

3n2  n^  Sin  a/2  COS  1/2(0 m  +  a)  -  Sin  l/2(9m  ♦  a)  COS  a/2 
30  ’  2  .  sin2  a/2 

3n2  nj  sin  (a  +  0m)/2 
30  2  Sin2  «/2 

Therefore  the  total  derivative  is 


dn„ 


Sin  1/2(6  +  or) 

_  Jut _ 

sin  a/2 


dnj  + 


hj  COS  1/2(6 m  + 


or) 


2  sin  a/2 


d6  + 
m 


nlBln6m/2 

2  _,_2 


dor 


sin  or/2 


At  this  stage  it  is  necessary  to  make  some  estimates  of  6m  and  a  because  these  enter  into  each 
of  the  partial  derivatives.  This  is  done  by  assuming  nominal  values  (4  and  1,  respectively)  for 
the  indexes  of  germanium  and  air.  Then 

sin  1/2(6  +  a) 

4  =  sin  a  /2 - 


Because  this  equation  has  no  unique  solution,  some  judicious  choices  must  be  made.  U  6^  +  a 
is  180°,  then  a  must  be  about  30°  and  6_  is  150°.  This  is  usually  too  large  a  dispersion  angle 
for  most  experimental  arrangements;  the  prism  angle  would  be  sufficiently  large  to  make  the 
optical  path  sufficiently  thick,  so  that  there  is  too  much  absorption.  If  the  prism  angle  is  10°, 
then  the  minimum  deviation  angle  is  30°;  this  is  more  reasonable.  By  the  original  considera¬ 
tions,  one  has 

dn^  sin  1/2(6  m  +  a) 

IHj  *  sin  a/2  *  * 


dn2  nj  cos  l/2(6m  ♦  a)  x /0.9401 /0.470\ 
86  m  ~2  - - sin  a/2  "  2V0.087/  \OV087 / 


5.402 


?n2  tin  6^/2  ^  q.259  b 
2  sin2  a/2  2(0, OUT)2 


0.1295  0.1295 


(0.087) 


2  0.007569 


17.10 
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Since  is  not  strongly  dependent  on  temperature,  the  assumptions  permit  numerical 
evaluation  of  the  equations  written  above  which  yields 

dn,  =  4dn,  +  5.402 d6^  +  17.10  da  +  2.68  x  10'4dT 


The  angles  6  and  a  are  expressed  here  in  radians.  Therefore  6  must  be  measured  to  2  x 

10  mrad  and  a  must  be  measured  to  5  x  10  mrad  to  obtain  1  part  in  the  fifth  decimal.  In 

more  familiar  units,  6  must  be  measured  to  about  one  half  second  of  arc  and  a  measured  to 
m 

one  eighth  second  of  arc. 


The  question  of  wavelength  uncertainty  remains.  The  dispersion  curve  for  any  optical 
material  is  nonlinear,  (if  it  were  linear,  one  index  measurement  at  each  of  two  wavelengths 
would  suffice  for  an  spectral  measurements  and  this  research  problem  would  be  almost  trivial.) 
Therefore,  errors  in  wavelength  are  proportional  to  index  errors,  but  the  "constant"  of  pro¬ 
portionality  is  a  function  of  wavelength.  In  the  region  around  2  p  (where  dn/dX  is  large),  errors 
in  wavelength  are  serious;  at  larger  wavelengths  they  are  not.  By  taking  published  values,  one 
can  iearn  something  about  the  permissible  values.  From  the  data  [5],  values  of  range  from 
0.0001  at  about  15  p  to  0.1  at  about  2  p.  A  reasonable  midband  value  of  0.01  yields 


||  dX  »  0.01  dX 

Therefore  dX  must  be  less  than  0.001  in  this  region;  in  fact  a  curve  can  be  constructed.  Some 
of  the  points  are  listed  in  Table  II.  The  worst  case  requires  spectral  accuracy  to  2.5  A,  and  this 
Is  short  of  3  p  where  there  are  enough  emission  lines. 

3.1.2.  PRISM  DEVIATION  METHOD.  In  this  method  a  beam  of  light  is  refracted  through  the 
prism  at  any  deviation — not  necessarily  minimum  deviation.  The  geometry  of  the  prism  provides 
the  first  equation: 


8'j  +  S' 


*  8 


Snell's  law  provides  the  next  pair  of  equations: 

nj  sin  ■  nj  sin  0'j 

nj  sin  $2  «  ng  sin  0r2 

Now  a  series  of  substitutions  puts  these  equations  in  a  form  that  simplifies  solution  of  n^/iy 


Institute  of  Science  and  Technology 


fho  University  of  Michigan 


ftj  sin  ■  n2  sin  (a  *  0'2) 

■  n2  (sin  a  cos  0'2  -  cos  a  sin  0'2) 

■  n2  join  a  (l  -  six?  O' -  cos  a  sin  0’2J 


Sin  0'o  ■  —  Sin  0- 
1  "2  2 


Sin  $l  a  sin  a  jj^nj/nj)2  -  sin2  *  cos  a  Sin  02 

2  o 

,  ,  sin  0,  +  2  sin  0,  cos  a  sin  0„  +  cos  a  sin  0_ 

(n2/n.)z  a  sin2  0  +  -  1  *•  .  2  2 

21  2  s  in2  a 

2  2 

0  Sin  0-  +  sin  0,  +  2  sin  0,  cos  a  sin  0. 

(„  /n  )2  =  2 - 1 - ,  1  - — 2 

2  1  sin2 a 

n2  *  -^a  (sin2  0j  +  sin2  ©2  +  2  sin  0j  sin  02  cos  aj1^2 

This  is  the  "simple"  form  for  n^  (at  least  it  is  an  explicit  representation  of  hj  or  i^/n^  in 
terms  of  prism  and  refraction  angles).  The  total  derivative  of  (a  function  of  the  angles  0^ 
and  02),  a,  and  the  index  n^  must  now  be  found: 

3n,  3iu  ?n_  3iu 

**2  *1  +|0j^ +W2d°2  +  W*9 


The  evaluation  proceeds  as  follows: 


m  iih  (ita*  °1 +  ,in2  V  2  ita  h  w  “J1 


0j  cos  0j  4-  cos  0j  sin  Sj  “•  °) 


sin  a  ^sin2  0j  fin2  02  +  2  sin  0j  fin  02  cos 

3n2  n^(aln  0j  cos  02  +  fin  0j  cop  0j  cos  a) 

992  sin  a  ^sin2  0^  +  sin2  02  +  2  pin  0j  pin  (j  cos  oj1^2 
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Now  it  is  necessary  to  make  some  assumptions  about  the  actual  experimental  arrangements. 
Assume  that  n^  is  1,  n^  is  4,  and  9^  and  02  are  20.4°.  Numerical  substitution  shows  that 

^2 

an. 


an, 


2n„ 


3ft* 


*  4, 


jf  =  5.37,  =  5.37,  and  —  =  22.7.  Therefore, 

1  2 


dn* 


4  dnj  +  5.37  d01  +  5.37  dflj  +  5.37  de^  +  22.7  do 


The  results  are  that 

n,  must  be  measured  to  2.5  x  10  ; 

*  _ 

91  and  0_  must  be  measured  to  1,86  x  10  ; 

1  z  -7 

a  must  be  measured  to  4.4  x  10  . 

With  the  assumptions  made  above,  this  method  is  about  equivalent  to  minimum  deviation 
because  and  9^  were  assumed  equal;  this  is  not  always  required,  however.  In  fact,  the 
Situation  is  such  that  as  gets  larger,  so  does  02,  and  the  error  increases.  This  is  not  par¬ 
ticularly  important,  since  the  experiment  would  be  set  for  minimum  deviation,  but  no  care 
would  be  taken  to  determine  that  Such  was  the  case.  The  error  analysis  for  minimum  deviation 
IS  then  essentially  correct. 


3.1.3.  TWO- POSITION  PRISM  METHOD.  The  Servo  Corporation  of  America  reports  a 
scheme  for  measuring  refractive  index  that  incorporates  a  reflecting  side  on  the  prism  [6]. 

The  arrangement  is  Shown  in  Figure  2. 

The  prism  angle  is  a,  the  rotation  angle  is  e,  and  their  sum  is  0.  The  deviation  angle  is  6; 
therefore  the  angle  relationships  are: 

6  -  -0g  +  Sj  0  »  or  +  « 

V*a  vfl  j 

j 

The  Snell  relations  are  j 

i 

1 

nj  sin  0j  m  n2  sin  0'j  *  ] 

I 

n^  sin  «  n2  sin  0’3  1 

•  I 

Therefore  j 

nj  Bin  0  ■  n2  sin  0'3  t 

From  the  trigonometry  of  the  triangles,  j 

10  1 
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’s  ****** 

W* 


When  these  equations  are  combined,  the  following  useful  relation  is  Obtained: 

S'3  *  2a  * 


therefore 


n^  Sin  0  *  ng  Sin  (2a  *  9’^) 

I"  /"l  8ln  ®l\l 

*  ng  Sin  2a  *  arcsln  ^  n  M 
Simplification  is  accomplished  thus: 


h« 

sin  0  •  -p  (sin  2a  cos  0’j  -  coS  2a  sin  S’j) 


sin  9\  ■  —  sin  0, 
1  n2  1 


c.' » 

1  n2i\v 


nl/2 

»J 


r-  a? 


sin  0  -  sin  2a 


•  sin3  9. 1  -  cos  2a  sin  (0  ♦  0) 


f* 


Solving  for  hg/nj  is  accomplished  thus: 

sin30  +  2  sin  0  sin  (0  +  6)  cos3  a  *  cos3  2a  sin3  (0  *  #)  *  sin3  2a 

3  »  .2 _ _ ._2 


ft)-*1”’9] 

n3\  _  sin3  0  +  2  sin  0  sin  (0  ♦  0)  cos  2a  +  cos3  2a  sin3  (0  ♦  0)  ♦  sin*  2a  sin3  (0  »  0) 

sin  2a 

n2  sin3  0  -f  2  sin  0  sin  (0  +  0)  cos  2  a  »  sin3  (0*0) 


G5 


sin  2a 

An  error  analysis  is  then  made  for  this  technique  in  the  usual  way. 


11 


3n  b _ sin  0  cos  {0  +  6)  cos  2a  +  sin  (0  +  0)  coa  (0  + 

sin  2a  [sin2  0  +  2  sin  0  sin  05  +  6)  cos  2a  +  sin2  (0  +  6)] 

cos  (0  +  6)[sin  0  cos  2a  +  sin  (0  +  0)] 

- - 5 - «  . 1 

Sin  2a  [sin  0  +  2  sin  0  sin  05  +  6)  cos  2a  +  sin  (0  +  6)] 


Numerical  evaluation  is  accomplished  by  assuming  that 


a  *  159 
0  *  3Q9 


«-  59 
2a  -30° 


fj  ■  3.029 

fi-™ 


Because  Oj  is  a  function  of  temperature  (besides  being  a  function  of  the  geometry),  and 
because  n^  is  not  strongly  temperature-  or  pressure-dependent,  the  result  is 

i>2  *  wij  +  f(T) 
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where  T  is  temperature,  then 


* 


dn2 


3n. 


dn 


an.  an 

I  "aiT  dn  +  iTr  ** 


ans 

ndn^  +  frjdn  '*  -pp 


dT 


^1*2  >1 

For  germanium  and  air,  the  index  values  are  n  *  4,  n^  *  1,  and  -pp  «  2.68  X  10  . 


Therefore 


dn2  =  4  dnj  +  3.025  d$  +  7.93  da  +  1.268  d6  +  2.68  X  10*4dT 

-6 

As  before,  n.  must  be  measured  to  an  accuracy  of  2.5  x  io  and  t  must  be  measured  to  o.037°C; 
*  -8 

6  must  be  measured  to  3.3  xio  rad  (Or  about  6.68  second  of  are);  a  must  be  measured  to 

-6  -6 
1.26  x  10  rad  (or  about  6.27  second  of  arc);  6  must  be  measured  to  7.89  x  io  rad  (or  about 

1 .63  seconds  of  arc). 


3.2.  REFLECTION  METHODS 


In  the  early  1950's  it  was  found  that  the  optical  constants  of  materials  with  metallic  proper* 
ties  could  be  determined  by  measuring  the  reflection  coefficient  and  the  phase  relationship 
between  the  two  components  of  the  reflected  electric  vector — one  component  parallel  and  the 
other  perpendicular  to  the  plane  of  incidence  (see  Figure  3).  There  are  basically  three  ways 
in  which  data  obtained  in  this  way  can  be  used  to  determine  the  optical  constants:  the  Ditchburn 
method,  the  Tousey  method,  and  the  Robinson  method.  These  methods  and  their  limitations  are 
described  below. 


f. 


3.2.1.  ditchburn  [7]  METHOD.  If  the  reflecting  surface  is  free  of  any  contamination,  one 
one  can  relate  the  phase  difference,  A,  between  the  two  different  orientations  of  the  electric 
vector  and  the  ratio  of  their  amplitudes  Rp/B„  (parallel  to  perpendicular)  to  the  optical  constants 
of  the  material.  This  relationship  is 

1  -  e^tan*  ,/n2  -  k2  -  sin2  0  -  tank 

1*  extant  *  W#Wd 


where  tan2  *  «  e  is  the  angle  of  incidence,  and  the  complex  index  of  refraction  is  tl  >  n  -  lk. 


Separating  the  real  and  imaginary  parts,  one  obtains 
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n2  *  k2  » tan2  a  Bln2  4 


COS2  24 
L  (1  * 


2*  -  sin2  2*  Bln' 


aln  24 f  cos  A) 


—  +  Sin2  a 


► 


and 


2  nk  » tan2  a  Bln2  a 1 


slA  4*  sin  A 


(1  +  sin  2*  cob  A)' 


2 


Therefore  n  and  k  can  be  calculated  from  these  last  two  equations,  on  the  basis  of  experimental 
values  of  A  and  4'. 

The  phase  A  depends  upon  0  and  k:  a  is  small  for  small  f,  and  gradually  increases  from 
0 o  to  iso®  as  a  changes  from  o  to  00°.  As  k  decreases,  the  region  of  transition  of  A  from  low 
to  high  values  becomes  more  pronounced;  when  k  is  aero,  A  is  a  step  function;  i.e.,  o°  to  180° 
at  the  Brewster  angle.  Thus  when  k  is  low,  the  phase  angle  is  difficult  to  determine  (possibly 
only  to  ±o.l°);  the  accuracy  of  the  a  determination  from  geometrical  considerations  is  about 
1  arcsec.  The  error  analysis  is  accomplished  in  the  usual  way.  First,  the  two  equations  above 
are  solved  for  n: 

4  2  Itan2  4  Bln2  4  (cos2  2*  -  sin2  2*  sin2  A)  .  2  71  tan4  4  Bln4  4  sin2  4»  sta2  A 

“  '"  I  J  ««. .»»«»*)' 

The  following  assumptions  are  made:  4  =  60°,  ♦  “  26.5°,  A  *  6°,  n  m  4. 

The  value  for  3n/d4  is  then  found  to  be  0.1.  The  choice  of  parameter  values  is  based  upon 
a  reasonable  geometric  arrangement  and  the  theory  presented  by  Dttchburn.  The  variation  in 
n  due  to  changes  in  4>  can  be  found  in  the  same  way:  0n/34>  is  -0.118.  Then  3n/9A  is  found  to 
be  0.0338.  Thus 

dn  -  0.1  d4  +  0.034  dA  -  0418  d* 


H  4  can  be  measured  to  one  arcsec,  A  within  one  urcmln  and  ♦  to  0.03  urcmln,  then 

dn  *  0.0011178 


This  implies  that,  for  wavelengths  shorter  than  1.5  p,  the  method  permits  determining  the 
constants  with  an  accuracy  of  10*2  if  the  effects  of  surface  films  such  as  oxides  and 
surface  irregularities  are  negligible.  However,  the  work  of  Avery  [8]  and  rise  that  of  Archer 
[8]  indicates  that  the  best  accuracy  obtainable  at  these  shorter  wavelengths  is  about  2%.  It  is 
therefore  unrealistic  to  assume  that  the  method  will  be  applicable  to  the  longer  wavelengths 
with  sufficient  accuracy. 


* 
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In  order  for  dn  to  be  -10  ,  the  accuracy  desired,  d<f>  must  be  no  greater  than  10  degrees, 

a5  .5 

d*  must  be  no  greater  than  10  degrees,  and  dA  must  not  exceed  3  x  10  degrees.  In  terms 
of  arcseconds,  the  accuracy  must  be  at  least 

”d <t>  <  0.36  arcsec 
Kd*  <  0.36  arcsec 
dA  <  1.08  arcsec 

Having  to  know  dd  this  accurately  could  be  quite  a  problem,  since  surface  irregularities 
could  cause  errors  that  exceed  the  required  precision.  Also,  the  Incident  beam  must  be  col¬ 
limated  very  carefully.  Polishing  the  sample  strains  the  surface;  this  then  necessitates  careful 
annealing  techniques.  These  surface  strains  have  been  known  to  introduce  inaccuracies  of  as 
much  as  5%  in  refractive  index  measurements  of  PbS  by  reflection  techniques.  To  date,  the 
error  in  determining  4  (due  to  unavoidable  experimental  problems)  amounts  to  about  3.6  arcsec, 
or  an  order  of  magnitude  greater  than  the  minimum  acceptable  value  of  0.36  arcsec. 

The  values  given  above  are  for  \  <  1.5  p  or,  in  a  region  of'  relatively  large  extinction  co- 

-7 

efficient,  k  >  0.5.  The  extinction  coefficient  in  the  region  of  Interest  is  around  10  .  This 
renders  the  determination  of  n  a  difficult  task  at  best. 


Another  problem  encountered  when  using  this  method  to  determine  refractive  index  is  the 
effect  of  film  formation  on  the  surface  of  the  sample.  Coatings  such  as  oxides  and  contaminants 
are  commonplace  and  extremely  difficult  to  avoid.  According  to  Archer  [9],  both  A  and  d  are 
affected  by  the  presence  of  surface  films.  In  his  measurements  a  10- A  oxide  film  on  Ge  caused 
a  2%  to  30%  error  inkanda0.6%to  5%  error  in  n.  Thin  oxide  films  obviously  can  greatly  disturb 
the  measurements  in  the  second  significant  figure  of  the  refractive  index.  This  is  quite  reas¬ 
onable,  since  the  measured  refractive  index  in  this  method  is  simply  the  index  of  the  sample 
surface,  inasmuch  as  the  reflection  mode  is  being  used. 

3.2.2.  TOUSEY  [10]  METHOD.  This  method,  too,  employs  reflection  measurements  made 
for  two  different  orientations  of  the  electric  vector. 

R  2  2  2  2 

p  _  a  +  b  -  2a  sin  d  tat,  ft  +  sin  d  tan  0 
-  .  ~  o  V 

n  a  +  b  +  2a  sin  d  tan  d  +  sin  d  tan  d 

where 

2a2  =  {[(n2  -  k2)  -  sin2  d]2  +  4n2k2}1/2  *  [(n2  -  k2)  -  sin2  d] 
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2b2  =  {[(n2  -  k2)  -  sin*  4>]‘  +  4n2k*}1/2  -  ((ft*  -  k2)  -  Bin2  0] 


2  a12  .  aJM2 


Now  the  measurement  of  ^  at  two  angles  of  incidence  will  provede  two  equations  with  two 

n 


unknowns: 


Rj-fjfn.k) 

R2  *  f2  (h,  k) 


These  equations  are  solved  graphically.  Tousey  states  that  the  uncertainty  in  n  is  about 
0.01.  This  is  10  times  worse  than  the  maximum  allowable  uncertainly. 


3.2.3.  ROBINSON  [11]  METHOD.  Just  as  the  phase  and  attenuation  of  the  Output  of  an 
electrical  network  are  interdependent,  so  is  the  phase  and  attenuation  of  the  reflected  electric 
vector.  Thus  one  can  use  reflection  at  normal  incidence.  The  Fresnel  equation  for  the  ampli¬ 
tude  of  reflected  light  at  normal  incidence  from  a  flat  surface  is 


The  reflection  coefficient  R  is 


,  (n  -  ik)  -  1 
(n  -  ik)  +  1 


Irle 


id 


R  =  Irr 


If  Irl  is  known  over  the  entire  wavelength  spectrum,  then  0  at  any  signal  frequency  can  be 
determined  from  the  Kramers-Kronig  relation 


The  success  of  the  method  depends  on  the  fact  that  negligible  error  results  from  a  lack  of 


knowledge  of  all  parts  of  the  frequency  spectrum.  This  technique  has  been  used  satisfactorily 
for  the  wavelength  interval  where  the  extinction  coefficient  is  quite  high  (i.e.,  on  the  short  wave- 

_7 

length  side  of  the  absorption  edge).  In  the  region  of  interest,  however  k  is  10  and  the  im¬ 
aginary  part  of  the  Fresnel  coefficients  vanishes;  the  result  is  that  the  refractive  index  meas¬ 
urement  reverts  to  a  determination  of  n  by  using  the  Brewster  angle  00: 


n  *  tan  0B 


Again  the  problems  of  surface  films  and  incident  angle  determinatives  are  encountered. 


3.3.  INTERFEROMETRIC  METHODS 

j  The  essential  idea  of  interferometers  involves  dividing  the  light  from  a  source  into  two  or 
'more  beams  which  are  teen  superposed;  tee  irradtance  in  tee  region  of  superposition  varies  from 
point  to  point  between  maxima  (which  exceed  tee  sum  of  tee  irradlances  in  tee  beams),  and 
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minima  (which  may  be  zero).  Several  different  Instruments  are  available  which  will  divide  a 
single  beam  in  order  to  exhibit  the  interference  pattern  [12] .  The  mathematics  of  the  inter¬ 
ference  phenomena  is  not  dependent  on  any  particular  instrument;  thus  the  conditions  for  fringes 
to  be  formed  will  be  discussed  in  a  general  way. 

The  irradiance  has  been  defined  as  the  time  average  of  the  amount  of  energy  which  crosses, 
in  unit  time,  a  unit  area  perpendicular  to  the  direction  of  the  energy  flow.  Then  the  total  ir¬ 
radiance  for  two  beams  propagating  in  the  same  direction  is 

H  -  Hj  +  H2  +  SyH"lt“  cos  6 
where  0  is  the  phase  difference  between  the  two  waves. 


Evidently  there  Will  be  maxima  of  irradiance  where 

lei  -  0,  2v,  4 1 . 

and  minima  of  irradinace  where 

lei  M.ji,..., 

V  a  plane  parallel  plate  of  transparent  material  is  illuminated  by  an  infinitely  distant  point 
source  or  a  collimated  beam  (see  Figure  4),  some  point  P  on  the  same  side  as  S  is  reached  by 
two  rays — one  reflected  at  the  upper  surface  and  the  other  at  the  lower  surface — so  that  there 
is  an  interference  pattern.  The  two  rays  from  S  to  P  travel  different  optical  paths  having  a 
difference  of  AL.  Then  the  corresponding  difference  in  phase  is 

±6  *  ^  AL  ±  n 

where  ±w  is  the  phase  advance  or  retardation  from  reflection.  From  the  geometry  in  Figure  4, 

AL  >  n'(AB  +  BC)  -  nAD 

where  n*  and  n  are  the  refractive  indexes  of  the  plate  and  the  surrounding  medium.  H  t  is  the 
r  thickness  of  the  plate,  and  I,  V  are  the  angles  of  incidence  and  refraction,  then 

AB  *  BC  « t/cos  P 

»  AD  ■  AC  sin  I  *  2t  tan  P  sin  P 

n  sin  I  m  n'  sin  P 


then 


AL 


2n't 


-  sin2P 
cos  P 
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Therefore 


There  are  bright  fringes  when 


There  are  dark  fringes  when 


AL  =  2n't  cos  I' 


±6  =  ^  n't  cos  I'  ±  it 


±6  =  ±2mir 


±C  =  ±(2m  +  i)# 


This  corresponds  to  bright  fringes  when 

2n't  Cos  I'  *  {m  +  “Jx 


There  are  dark  fringes  when 


In  each  case 


2n’t  cos  T  =  mX 


m  =  0,  1,  2,  .  .  . 


A  given  fringe  is  characterized  by  a  constant  value  of  I’  and  therefore  a  constant  value  of  I 
because 


n  sin  I  -  n'  sin  I' 


The  condition  for  circular  fringes  is 


cos  !' 


l/Hl 


sin'  ! 


3.3.1.  SIMPLE  INTERFEROMETRIC  METHOD.  The  refractive  index  is  found  (onthe basis 
of  the  theory  above)  to  be 


*  "2 


where 


m  c  m  or  mix 
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The  differential  error  may  be  written  as 


In  performing  the  actual  experiment,  the  angle  of  incidence  will  be  almost  zero  and  the  refractive 
index  of  air  is  nearly  one.  The  wavelength  will  be  allowed  to  vary  from  1  to  is  jx.  The  values 
Chosen  are  n'  =  4,  t  =  5  mm,  and  1  =  4  arcmin.  The  allowable  errors  are  then  dl  -  3  arcmin,  dt 
■  i.25  A,  dx  *  id"5  ^  dmQ  *  o.oi,  and  do  *  id. 

One  way  to  evaluate  the  usefulness  of  this  method  is  to  examine  the  restriction  placed  on  the 

parameters.  The  severest  restriction  placed  on  this  method  is  the  measurement  of  length, 

(i.e.,  parallelism).  Obtaining  a  narrow  source  is  also  a  problem.  The  use  of  a  line  source  (e.g., 

a  hydrogen  gas  discharge)  is  ruled  Out  since  there  are  no  sources  of  this  type  in  the  spectral 

region  of  interest,  if  there  were  an  ideal  source,  still  another  problem  would  be  introduced: 

4  3 

determining  m  (4  x  Id  <  m  <  2.5  x  10  ),  the  order  of  interference.  If  the  desired  resolution 
can  be  obtained  and  the  number  Of  fringes  Can  be  determined,  the  next  step  would  be  to  determine 
the  thickness  to  *125  A;  this  IS  impossible,  as  was  shown  in  the  section  on  reflection.  Also, 
the  two  faces  of  the  sample  must  not  make  an  angle  greater  than  0.001  second  of  arc;  they  must 
be  extremely  parallel,  because  of  the  restriction  on  thickness. 


The  refractive  index  of  air  presents  no  great  problem. 

The  next  error  is  that  concerning  the  angle  of  incidence,  which,  like  the  refractive  index 
of  air,  does  not  introduce  any  great  problem.  In  fact,  at  small  angles  the  angular  resolution 
required  decreases  and  would  be  negligible  for  an  angle  of  (for  example)  less  than  1°  off  normal. 
This  happens  because  the  sine  function  decreases  at  small  angles  and  is  the  major  contributor 
to  the  angular  error. 

3.3.2.  ROTATING  SAMPLE  METHOD.  Another  experiment  could  be  constructed  which 
employs  two  plates;  one  plate  rotates  with  respect  to  the  other.  The  refractive  index  of  a  trans¬ 
parent  plate  and  its  dispersion  can  be  determined  by  the  use  of  white  and  monochromatic  light. 
The  determination  of  the  dispersion  is  based  upon  the  shift  between  the  true  and  the  apparent 
position  of  the  center  of  the  system  of  fringes  formed  by  white  light  [13], 

A  plate  is  cut  in  half,  and  one  of  the  halves  is  rotated  slowly  about  an  aids  perpendicular 
to  the  beam;  the  fringes  begin  to  move  because  of  the  increase  in  the  optical  pathlength  (see 
Figure  5),  From  the  geometry,  the  difference  in  optical  path  can  be  written  as 
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these  expressions  assume  that  the  instrument  is  in  perfect  alignment,  that  the  plates  are 
parallel,  and  that  the  source  is  Collimated.  Obviously  these  conditions  are  not  attainable  in 
any  real  experiment.  The  limitation  on  the  accuracies  are  discussed  in  Section  3.3.1. 

By  setting  n'  =  4,  the  following  resultB  Can  be  obtained:  I  a  3°.  m  -  10  at  X  *  2  a.  and  t  * 

©  A  -B 

5  mm.  then  dl  ■  0,5  arcsSc,  dt  *  100  A,  dX  *  2  x  10  a,  dm^  *  2  x  10  fringe,  and  th»a4r  * 

5  x  10’7. 

the  limitations  by  this  method  are  the  same  as  in  Section  3.3.1  except  that  the  number  of 
fringes  to  be  counted  is  on  the  order  of  10;  the  method  still  does  not  seem  physically  realizable. 


3.3.3.  RELATIVE  INDEX  MEASUREMENT.  Another  method  can  be  devised  to  reduce  the 
number  of  fringes  that  must  be  counted.  If  the  refractive  index  is  determined  at  one  wavelength, 
it  can  easily  be  determined  for  any  other  wavelength  by  use  of  the  difference  relation 

n'(Xj)  -  n'(Xjj)  *  jt^os  V  (X1  ”  V 

The  position  of  normal  incidence  (cos  T  =  1)  could  be  found  by  rotating  the  plate  until  a  null 
occurs  (l.e.,  when  the  fringe  pattern  starts  to  expand  for  either  direction  of  rotation).  At  this 
point  the  other  parameters  could  easily  be  determined,  except  for  n'(Xj). 

The  errors  can  be  determined  by  the  differential  method: 

(x.  - x  ) 

dn>(x?).dn-(xJ).|j^rfd(Am) 
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Am 
2t  cos 


P^VV1 


AmfXj  -  \2) 
2t  cos  I' 


dt 

t 


AmfXj  -  X2) 
2t  Cos  V 


tan  I’  dl 


To  determine  the  total  error,  the  following  data  will  be  used:  Xj  *  Xg  e  1/2  p,  t  *  5  min,  and 

Am  =  2.  It  is  assumed  that  one  knows  h'(X. )  to  ±0.00001  and  that  the  beam  impinges  at  normal 

incidence  (cos  I*  I);  then  for  a  =  5  x  io  ,  t  =  2  x  10  m,  I  *  1.2  x  10  rad,  the  allowable 

m  _« 

uncertainties  are  d6m  =  0.01  fringe)  6 AX  =  10’  p,  and  dt  «  0.01  mm. 

This  method  (relative  index  measurement)  seems  to  be  the  best  of  those  considered,  but 
depends  on  an  initial  determination  of  the  refractive  index.  The  physical  parameters  are  well 
within  experimental  procedures.  The  variation  of  the  index  could  easily  be  checked  for  the 
wavelength  region,  one  fringe  at  a  time,  and  a  complete  set  of  data  could  be  obtained.  The 
temperature  variation  could  also  be  found  if  the  sample  were  heated  while  the  fringe  pattern  was 
observed.  This  method  would  also  Offer  a  check  on  the  variation  of  the  index  measured  by 
another  method. 

From  the  experimental  standpoint,  two  other  problems  are  introduced:  alignment  of  the 
instrument,  and  heat  transfer  in  the  Sample.  There  are  straightforward  procedures  available 
for  the  alignment  of  interferometers,  but  these  are  delicate  and  time-consuming. 

If  the  absolute  index  can  be  determined  by  some  other  method,  the  interferometric  technique 
is  very  good  for  determining  any  variation  in  the  index  as  a  function  of  wavelength  and/or 
temperature. 


3.4.  SPHERE  METHOD 

This  method  of  measurement  utilizes  multiple  internal  reflections  through  a  solid  Ge  sphere. 

A  light  ray  incident  upon  a  homogeneous,  isotropic,  transparent  sphere  with  an  index  n  is 
refracted  into  the  sphere  and  undergoes  p  - 1  internal  reflections  with  a  portion  of  t!.e  ray  being 
refracted  outside  at  each  internal  reflection.  The  total  deviation  of  the  ray  from  its  initial 
direction,  $,  is  given  by 

$  =  21  -  2 pi'  +  (p  -  l)» 

where  I  is  the  angle  of  incidence  and !'  is  the  angle  to  the  normal  after  refraction  (obtained  from 
Snell's  law).  A  particular  ray  which  satisfies  the  equation 
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is  called  a  "limiting  ray."  This  equation  can  be  satisfied  if  p  £  2  and  n  <  p.  It  can  be  shown 
that  *  is  a  minimum  for  a  limiting  ray;  when  the  emerging  light  is  viewed,  therefore,  a  concen¬ 
tration  of  light  is  observed  close  to  the  limiting  ray  and,  on  one  side  of  the  bright  band,  there 
will  be  darkness.  (The  bright  side  IS  actually  a  diffraction  pattern.) 

Calculation  also  shows  that  the  direction  of  the  limiting  ray  depends  on  the  refractive 
index  of  the  sphere  through  the  relation 

d$  2  /p^  - 

Thus  the  determination  of  0  allows  calculation  of  the  refractive  index. 

The  advantage  of  this  method  is  that  the  change  in  *  for  a  given  change  in  n  IS  much  larger 
than  is  the  change  in  the  minimum  deviation  angle  of  a  prism  for  the  same  change  in  n. 

In  order  to  design  an  experiment  to  measure  n  to  within  ±0.00001,  the  required  angle¬ 
measuring  accuracy  and  the  Size  of  the  sphere  required  must  be  determined.  A  third  consider¬ 
ation  is  the  amount  of  energy  arriving  at  the  detector  for  a  given  amount  incident  upon  the  sphere. 

A  detailed  examination  of  the  experimental  procedure  is  necessary  to  establish  the  accuracy 
requirement  on  angle  measurements.  The  position  of  a  limiting  ray  is  determined  by  measuring 
the  positions  of  two  minima  in  the  diffraction  pattern  associated  with  the  limiting  ray.  The 
angular  difference  in  seconds  in  the  two  minima  is  equal  to  a  constant  P  multiplied  by  the 
difference  in  the  v  values  for  the  two  minima  given  in  Table  IV  of  Walther's  thesis  [14] .  P 
can  be  determined  and  used  with  one  of  the  v  values  and  the  corresponding  minima  position  in 
order  to  calculate  the  position  of  the  limiting  ray.  If  the  accuracy  in  measuring  the  position  of 
the  minima  is  a,  then  the  accuracy  of  the  position  of  the  limiting  ray  is  V'SI. 

It  is  not  desirable  to  have  to  measure  the  position  of  three  limiting  rays;  e.g.,  those  cor¬ 
responding  to  p  »  5,  8,  7,  From  these  data,  or  *g  *  can  be  obtained  and  n  calculated 

without  knowing  the  direction  of  the  incident  ray.  If  each  limiting  ray  is  known  to  have  an  accuracy 

of  V55,  then  the  sum  of  two  limiting  rays  will  be  known  to  KSa.  A  calculation  shows  that  A(*. 

-5  - 

*  *  )  must  be  4  arcsec  for  —  -  -  *  .  Therefore  a  -  <YE  •  1.6  arcsec.  Thus  angle  position 

measurements  must  be  accurate  to  1.6  arcsec.  If  the  limiting  ray  position  is  accurate  to  FSa, 

then 
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or 


Therefore 


PAv  -  A*  (VS  1.6), 


Ai>  « 


(1.7)(1.6)  1.7  x  1.6 

P  300 


*0.01 


A*  *  0.01  P 


An  » 


0.01  P 
d*/dn 


Walther  gives  P  as 


1/3,.,2/S 


ww 


(2.06266  X  10°) 


where 


h  ■ 


pV-1) 


and  the  numerical  factor  is  the  conversion  from  radians  to  arcseconds.  V  this  is  used  without 
the  numerical  factor  and  substituting  from  above, 


However,  for  ^  ■ -  X^-,  r  >  11.2  cm. 

Thus  in  order  to  measure  n  to  an  accuracy  of  two  parts  in  the  fifth  decimal  place  requires 
a  Ge  sphere  about  23  cm  in  diameter  and  the  capability  of  measuring  angles  to  1.6  seconds  of 
arc.  This  required  diameter  is  impractical  and  inconsistent  with  P  «  300  which  assumes  r  *  4 
cm. 

H  the  accuracy  requirement  is  relaxed  to-p- 1  x  io"B  and  the  measuring  accuracy  is 
retained  at  a  ■  1.6  arcsec,  then  the  radius  required  is  4  cm;  this  is  practical  and  consistent 
with  the  value  of  P  used.  Note  that  it  is  not  allowable  to  relax  the  measuring  accuracy,  as  this 
will  change  the  numerical  factor  and  again  an  lmpractically  large  sphere  will  be  necessary. 
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The  third  consideration  is  the  required  energy.  The  absorption  coefficient  of  Ge  IS  approxi¬ 
mately  1  cm"1  at  15  4.  The  pathlength  In  the  Ge  1S»  2(p  x  r).  Therefore,  for  p  ■  5  and  r  ■  4 
cm,  the  attenuation  factor  is  about  e"*°  *  4  x  10’1®.  Even  if  all  other  attenuation  factors  are 
neglected,  the  experiment  is  impossible  because  of  internal  absorption  within  the  Ge.  if  the 
sphere  size  is  reduced  to  the  point  where  a  signal  is  detectable,  then  the  accuracy  is  less  than 
l  part  in  the  fourth  decimal  place. 


4 

TRANSMISSION  MEASUREMENTS 

A  sample  of  germanium  which  was  submitted  to  the  weather  Bureau  by  a  vendor  exhibited 
very  poor  transmission.  The  sample  was  investigated  in  order  to  learn  why  and  to  become 
further  acquainted  with  the  engineering  problems.  Transmissions  of  sections  of  this  sample 
were  measured  and  the  mechanical  yield  properties  and  their  relation  to  the  transmission  in¬ 
vestigated. 


4.1.  STANDARD  TRANSMISSION  MEASUREMENTS 

After  the  Very  low  transmission  of  the  1-inch  thick  sample  was  measured,  it  was  cut  in  half 
and  the  new  surfaces  were  ground  and  polished.  During  the  surface  preparation,  care  was 
exercised  to  minimize  strains  and  surface  stresses  introduced  during  the  cutting  and  polishing 
operations*  The  final  surface  finish  was  microscopically  comparable  to  that  observed  on  the 
original  (1  by  0.75  inch)  faces  of  the  sample* 

These  new  sections  were  arbitrarily  designated  A  and  B,  and  the  transmission  characteristics 
of  each  were  again  determined.  Marked  differences  were  observed  in  the  transmission  peaks 
in  the  1.8-4  region.  Sample  A  showed  a  peak  transmission  of  about  13%,  whereas  the  peak  of 
sample  B  was  24%.  These  results  indicate  pronounced  structural  differences  through  the  cross 
section  of  the  original  sample,  and  different  numbers  of  free  carriers, 

These  variations  were  studied  further*  Sample  A  and  sample  B  were  each  cut  in  half, 

-  providing  four  sections  0.203  ±  0.001  inch  thick  after  polishing,  Figure  e  illustrates  the  trans¬ 
missions  of  these  sections,  and  relates  each  section  to  its  position  in  the  original  sample*  These 
results  show  that  there  is  an  increase  in  the  free-carrier  absorption  from  sample  B  to  A;  this 
indicates  that  this  crystal  was  grown  under  nonuniform  conditions.  Sample  B  appears  to  con¬ 
tain  fewer  impurities  than  the  other  sections.  Even  here,  however,  grain  boundaries  were  ob¬ 
served  on  the  edges  of  the  section  and  internal  strains,  lattice  dislocations,  and  chemical 
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impurities  could  account  for  the  decreasing  absorption  following  the  2- p  peak.  It  is  apparent, 
therefore,  that  this  polycrystalline  sample  contains  an  impurity  gradation  from  one  side  of  the 
window  to  the  other. 


4.2.  STRUCTURAL  YIELD  PROPERTIES 

These  properties  were  investigated  by  the  dynamic  spherical -indenter  technique  (described 
in  Reference  15),  which  has  been  applied  to  both  amorphous  and  crystalline  structures;  the 
characteristics  of  the  defects  produced  by  the  indenter  trace  line  are  related  to  the  internal 
energy  and  bonding  characteristics  of  the  material.  The  defects  in  glasses  so  produced  are 
sensitive  to  variations  in  ionic  substitutions  within  the  network  [16] . 

By  examining  the  detailed  variations  in  the  flaw  characteristics  as  a  function  of  applied  load, 
one  can  readily  distinguish  amorphous  from  crystalline  structures,  and  determine  the  relative 
response  of  the  given  structure  to  applied  stress.  The  mean  length  of  the  gross  slip  lines 
created  by  the  indenter  is  designated  as  and  is  taken  as  a  measure  of  the  relative  bond 
strength.  The  flaw  number,  designated  Fn,  is  a  function  of  the  rigidity  of  the  network.  Hie 
product  of  these  two  parameters  is  an  indication  of  the  overall  response  of  the  given  system  to 
applied  stress  and  is  designated  by 

V„*“b  «> 

where  Nb  is  a  relative  measure  of  the  number  of  bonds  disrupted.  The  relationship  between 
stress  and  the  radial  distance  from  the  point  of  loading  on  this  spherical  indenter  is  given  by 


where  P  is  the  applied  load  and  v  is  Poisson’s  ratio.  The  critical  stress  o  occurs  at  the  point 
where  those  flaws  which  were  initiated  at  the  trace  line  terminate;  e  is  found  by  substituting 

v 

F  i  for  r  in  Equation  2. 

ec  -  kP/Ff2  (3) 


where 


k  •  (1  -  2v)/2t 


(4) 


If  one  assumes  that  a  is  a  constant  for  each  type  of  structure  and  that  crystal  orientation  is 

~  . .  "  “  *  c  . 


given,  the  flaw  length  is  found  to  be  related  to  the  applied  load  by 

F|  -  (Voc)1/2P,/2 


(S) 
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Equation  5  predicts  a  linear  relation  between  and  the  square  root  of  the  applied  load  on  the 
indenter.  Thus  under  these  given  experimental  conditions,  the  rate  of  flaw  change  with  load  is 


p"l/2 


(6) 


From  Equation  6,  the  energy  of  bond  displacement  expressed  as  fiaw  formation  may  be  equated 
with  the  force  exerted  on  the  test  surface: 


(7) 


After  integration,  the  functional  relationship  between  flaw  length  and  load  is  found  to  be 


+  K 


(8) 


where  k  is  a  constant  and  is  physically  related  to  the  length  of  flaw  created  by  a  threshold 

stress  which  must  be  introduced  into  a  given  network  to  initiate  the  flaw  formation.  Equation 

8  has  been  experimentally  confirmed  [2]  and  utilized  in  the  analysis  of  germanium  crystals. 

From  Equation  8  and  the  empirical  relationships,  0  is  given  by 

c 

ff,  =  k/a2  (9) 


where  a  is  the  slope  determined  from  the  experimental  curve. 

The  results  of  applying  the  indenter  to  two  sections  Of  the  germanium  sample  are  shown 
in  Figure  7.  Each  point  on  these  curves  is  the  average  of  60  measurements  of  the  length  of  the 
gross  slip  or  defect  lines  formed  along  the  indenter  traces  on  the  polished  germanium  surface. 
These  effects  were  displayed  by  etching  the  sample  for  two  hours  in  a  solution  designated 
"number  2  etch"  [17] ,  Although  there  is  considerable  scatter  in  the  points  (see  Figure  7),  the 
results  generally  follow  the  linear  relationship  given  by  Equation  8.  The  mean  values  for  both 
the  flaw  lengths  (F()  and  the  flaw  number  (Fn)  are  given  in  Table  n  for  sections  B  and  C. 

The  critical  stress  can  be  calculated  from  the  curves  in  Figure  7  by  using  Equation  9. 

These  values  are  listed  in  Table  in.  The  values  of  the  calculated  critical  stress  show  a  general 
agreement  with  the  optical  transmission  characteristics  shown  in  Figure  6.  Section  B,  which 
shows  the  highest  transmission,  would  also  contain  the  lowest  impurity  level  Mid  fewer  structural 
defects  (e.g.,  slip  lines  and  dislocations).  Under  an  applied  stress,  section  B  would  permit  a 
greater  degree  of  localised  yield;  a  higher  stress  would  therefore  have  to  be  applied  before  the 
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gross  defect  lines  were  formed.  Section  C,  however,  evinces  a  very  high  impurity  level.  This 
would  make  dislocation  pinning  more  likely  and  would  inhibit  glide  mechanisms;  section  C  would 
therefore  exhibit  more  brittle  tendencies,  and  the  ac  value  would  be  predicted  to  be  lower. 
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FIGURE  4.  GEOMETRY  OF  INTERFERENCE  AT  A  PLANE  PARALLEL  PLATE 


FIGURE  5.  GEOMETRY  FOR  ROTATING  PLATE  METHOD 
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TABLE  n. 

WAVELENGTH  ERRORS 

X 

dn/dX 

(pH 

dX 

M 

(#) 

3 

0.04 

0.25  x  10‘3 

10 

0.0005 

0.20  X  10'1 

13 

0.000125 

0.80  x  10*1 

TABLE  m.  CALCULATION  OF  CRITICAL  STRESS 
FOR  GROSS  DEFECT  FORMATION 

Section  (pal) 

B  159,000 

C  48,000 


TABLE  IV.  INDENTER  FLAW  PARAMETERS  IN  SECTIONS  FROM 
A  GERMANIUM  CRYSTAL 
Section  B 


(grams) 

F|  (mm) 

F„  (number  per  mm) 

*b 

2600 

0.0325 

7.47 

0.243 

1600 

0.0276 

6.10 

0.170 

1000 

0.0220 

8.07 

0.177 

750 

0.0155 

7.14 

0.111 

500 

0.0133 

6.21 

0.083 

200 

0.0111 

4.43 

0.049 

Section  C 

(grams) 

Fj  (mm) 

Fn  (number  per  mm) 

Np 

2600 

0,0540 

10.87 

0.587 

1600 

0.0263 

8.62 

0.227 

1000 

0.0194 

6.49 

0.126 

750 

0.0169 

6.02 

0.102 

500 

0.0116 

5.81 

0.067 
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